Quantum evolution of the Universe from r = in the 
constrained quasi-Heisenberg picture 



S.L. Cherkasf and V.L. Kalashnikov$ 

f Institute of Nuclear Problems, Bobruiskaya 11, Minsk 220050, Belarus 

| Technische Universitat Wien, Gusshausstrasse 27/387, Vienna A-1040, Austria 

Abstract. The Heisenberg picture of the minisuperspace model is considered. The 
suggested quantization scheme interprets all the observables including the Universe 
scale factor as the (quasi) Heisenberg operators. The operators arise as a result of the 
re-quantization of the Heisenberg operators that is required to obtain the hermitian 
theory. It is shown that the DeWitt constraint H — on the physical states of the 
Universe does not prevent a time-evolution of the (quasi)Heisenberg operators and their 
mean values. Mean value of an observable, which is singular in a classical theory, is 
also singular in a quantum case. The (quasi) Heisenberg operator equations are solved 
in an analytical form in a first order on the interaction constant for the quadratic 
inflationary potential. Operator solutions are used to evaluate the observables mean 
values and dispersions. A late stage of the inflation is considered numerically in the 
framework of the Wigner-Weyl phase-space formalism. It is found that the dispersions 
of the observables do not vanish at the inflation end. 
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1. Introduction 

A focus of the article is the construction of the quantization scheme in which the 
observables (including Universe scale factor) are the time dependent operators. This 
allows to evaluate their mean values and dispersions. 

A variety of the quantization schemes for the minisuperspace model can be roughly 
divided in two classes: imposing the constraint "before quantization" [Tj and "after 
quantization" In the former the constraints are used to exclude the " nonphysical" 
degrees of freedom. This allows then to construct the Hamiltonian acting in the reduced 
"physical" phase space. 

The last schemes prefer imposing the constraint "after quantization". This leads 
to the Wheeler-DeWitt equation on quantum states of the Universe El- Our 
scheme can be considered among the last class because we use the Wheeler-DeWitt 
equation. However we supplement this equation with the system of the quasi- Heisenberg 
operators acting in the space of the solutions of this equation. Quantization rules for 
these operators are defined consistently with the choice of the hyperplane used for 
normalization of the solutions of the Wheeler-deWitt equation in the Klein-Gordon 
style. 

In the simplest case of an isotropic and uniform Universe filled with a scalar field, 
the minisuperspace equation contains only two variables: the scale factor of Universe 
and the amplitude of the scalar field. There is no an explicit "time" in the corresponding 
Wheeler-deWitt equation, whereas we are interested namely in the Universe evolution in 
time. This leads to various discussions about "time disappearance" and interpretation 
of the wave function of Universe jl]. Possible solutions like to introduce time along 
the quasi-classical trajectories, or subdivide Universe into classical and quantum parts 
were offered [5J. Such point of view can not be satisfactory. Ideally, time does exist 
independently regardless of whether we consider Universe quantum or classically. That 
forces to consider the quantum evolution of Universe without any need in foreign 
ingredients [HI El E] • Let us remind that this situation is analogous to that in string 
theory, where the constraint H = also exists. Nevertheless this constraint does not 
prevent an evolution of the Heisenberg operators X(cr, t) along r [TO] . 

In fact the constraint H = (in a theory of the constrained systems H is usually 
refereed as the super-Hamiltonian) tells nothing about whether the mean values of the 
Heisenberg operators evolve with time or not. This is defined by the normalization of 
the wave function. We can outline an issue in the following way. Let there is the wave 
function ip(a,x) dependent on two variables and obeying Hip(a,x) = 0. 

For evolution of some Heisenberg operator H we have 

< A(t) >=< i)\e iilT Ae~ ikr \i) > . (1) 

At first sight it seems that there is no evolution. However this is not the case. Really 
e~ lHT \ip >= \ip >, but we cannot write < ip\e tHT =< ip\. The point is that the wave 
function can not be normalized in the ordinary way / ip*(a, x)ip(a, x)dxda = 1 due to the 
constraint. In fact the function is unbounded along one of the variables. For instance, 
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let this is a variable. So if H contains differential operator like ^2, one cannot move 
A to the left by the habitual operation < i>\-^2 =< JjrV'l through the integration by 
parts. As a result, the assumption that < ip\H =< Hip\ = is wrong. 

Unfortunately this is an oversimplified picture. Reality was found to be more 
complicated. In the next section we show that the normalization of the wave function in 
the Klein-Gordon style permits evolution of the mean values of the Heisenberg operators. 
However, additional efforts are needed to obtain a genuine hermitian theory. 



2. Quantization of a particle-clock 

The problem of the quantum cosmology has many common features with those of the 
relativistic particle QjJ EH E2] • The action for the relativistic particle can be defined 
as: 



S = -m J y/-a%dT. (2) 

The equivalent form leading to Eq. (J2J) by means of variation of the lapse function 
e(r) is 



S = -j(e x x\- em 2 )dr, (3) 

where (— , +, +, +) signature is used. One more equivalent form resulting in (J3j) due to 
varying looks as 



(4) 



Using the re-parameterization invariance [TU1 ITT] we can choose the lapse function 
as e = 1/m. One can see from the last equation that the Hamiltonian is equal to 
H = M 2m and vanishes on the constraint surface H = after varying on e. The 
quantization procedure is based on the assumptions: 

r , x v \ = -igr, p = {i, p} = {i^, -< v>. 

Then the constraint becomes the Klein-Gordon equation. The commutators of the 
position and four-momentum operators with the Hamiltonian result in the Heisenberg 
equation of motion: 

= = f. = . (5) 

ar m ar 

Evident solution of the motion equation is 

X N T ) = x m + ^W T; jp( T )=pf l (0), 
m 

where x M (0) = x^ = {t,r} and p M (0) = Because before quantization we do not use 
any gauge fixing providing the time-dependent second class constraint fTJ, t(r) can be 
treated as an operator. We shall refer to t(r) as "physical time" and to r as "proper 
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time" . Whereas the physical time is operator, the proper time is some parameter "always 
running forward" like time in the Newtonian physics. 

One can check that the wave function satisfying the Klein-Gordon equation cannot 
be normalized through the integration over d 4 x j^j- However, it can be normalized 
through the time-like component of the conserved four-current: 



where e 2 (k) = k 2 + m 2 . To avoid an appearance of the states with the negative norm we 
have to choose only positive frequency solutions. For this wave packet the integration 
/ ip*(x)ip{x)d 4 x does not converge. Hence the wave function is unbounded along the 
t-variable. 



where "h.c." means the hermitian conjugate value. 
The definition © has the following properties: 

1) It is consistent with the normalization of the wave function if we choose A to be equal 
to the unit operator. 

2) It looks like as an expression for the mean value of the Heisenberg operator in the 
nonrelativistic picture when the operator acts on the wave function taken at the initial 
moment of time (note, that after integration over d 3 r we have to set t = 0). 

3) It has a natural property <A £} > =< A j^- >• 

4) It gives zero for the physical time, when the proper time is equal to zero. 

Let us emphasize the difference between our quantization scheme for the relativistic 
particle and e.g. the quantization schemes (in the ADM style with the some gauge fixing 
or in the Foldy-Wouthauthen representation of the Klein-Gordon equation) where the 
time is not an operator. Let consider some observer with an "exact" clock tracking 
a particle motion. The observer can measure, for instance, the particle coordinate 
and describe it by the operator r(t) to take into account the quantum fluctuations. 
There remain only technical problems like to find the Hamiltonian translating f(t) = 
e lHt r(0)e~ tHt etc. Let consider another situation, namely the observer has no a clock but 
there is an exact "clock" on the particle (e.g., this is the "radioactive" particle emitting 
photons every equal time pieces in a reference system connected with the particle). 
The observer having no own clock is forced to measure time by the particle clock. The 
measured time is the fluctuating value described by the operator t(r) = —r, because the 




Normalized wave packets have the form: 

f a(k)e- iei - h)t+ikr 




Let us define the mean value of some Heisenberg operator A(t) for the case of a 
free relativistic particle as 




(a) (b) 





Figure 1. (a) Observer describing a particle motion by own clock, (b) Observer 
describing a particle motion by " particle clock" . 



particle energy e is the fluctuating quantity for the wave packet. Thus our quantization 
scheme describes namely the particle clock (i.e. the particle supplied with the clock) 
rather than the "mute" particle. 

After averaging, the Heisenberg equations of motion result in 

T T 

< t(r) >=< e > — , < f (t) >=< r > + < p > — . (7) 
m m 

One can see that the physical time is proportional to the proper time. It is interesting to 

calculate the dispersion of the physical time for the Gaussian wave packet a(k) ~ e~ afc2 

with the - >> m 2 . An evaluation of the mean values of e = ij^ and its square gives: 

A _ f™e{k)e- ak2 k 2 dk ^ 2 „ 2 _ f™e 2 (k)e- ak2 k 2 dk _ 3 

<£> ~ ^e-^kUk <£ >_ f™e- ak2 k 2 dk ~2^' 

that results in 



< i 2 > - < i > 2 V< i 2 > - < i > 



2 




< t > < e> 

Thus, the "particle-clock" is a bad clock when the particle is localized in the spatial 
region less than the Compton wave length. Let us imagine what occurs when the 
particle possessing the electric charge is placed in the electric field. Let there is a cloud 
of the particles with some initial dispersion of the energy placed in the electric field. The 
particles begin to accelerate in the electric field and thereby they receive the energies 
much larger than their initial energies. As a result the relative energy dispersion becomes 
negligible. Similar picture holds for the quantum case so that the relative accuracy of 
the "particle-clock" increases. 

Let us remind that in fact we have used here the covariant proper time formalism, 
which goes down to Fock jT^] and Kramers This formalism is widely used for 

derivation of the Bargman-Michel-Telegdi equations for the particle spin motion in the 
external fields [15J. In this formalism we have to normalize the wave function through 
the Klein-Gordon scalar product instead of normalization through the integration over 
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d A x (on the other side this permits the time-evolution despite of Hip = 0). Absence 
of hermicity leads to the complex mean value of some observables. Although we are 
try to avoid a complexity by adding "/i.e." quantity in Eq. (JHJ) the negative dispersion 
appears when we evaluate, for instance, < t 4 (r) > etc. Thus we have seen that the 
imposing constraint to the state vectors is not sufficient to obtain a "good" theory in 
the framework of the covariant proper time formalism. In next section we shall propose 
solution for this problem. 



3. Quantum cosmology 

Let us start from the Einstein action of a gravity and the action for an one-component 
real scalar field: 

S = J d 4 x^R + J d A x^g[\{d^f - Vffl], (8) 

where R is the scalar curvature and V is the matter potential which includes a possible 
cosmological constant effectively. We restrict our consideration to the homogeneous and 
isotropic metric: 

ds 2 = N 2 (r)dr 2 -a 2 (r)da 2 . (9) 

Here the lapse function N represents the general time coordinate transformation 
freedom. For the restricted metric the total action becomes 



a 2 



s = J *< T W ( K - wiF) J + r aw, - aV W T > 

where K, is the signature of the spatial curvature. This action can be obtained from the 
following expression by varying on p a and p^ 

Varying on N gives the primary constraint 

*~ <«* 

After quantization [a,p a ] = —i, [<j),p^\ = % this constraint turns into the DeWitt 
equation Hi/j(a,4>) = 0. Looking at the constraint equation a desire may appear to 
modernize or remove it ^Hl • Apparently this implies (both on classical and on quantum 
levels) existence of some preferred system of reference. Although there are some 
logically consistent theories implying preferred system of reference, for instance, the 
Logunov relativistic theory of gravity [T7] giving an adequate description of the Universe 
expansion |18j . we shall keep to the General Relativity here and retain the constraint. 

Let us first consider the flat Universe {K = 0) with V(</>) = (corresponding 
Hamiltonian is H ). 

Procedure, which is invariant under a general coordinate transformation consists in 
postulating the quantum Hamiltonian |19j 

= ^g~*p„gH^Pvg~K (n) 
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where p^ = -ig \JL. g z = -i (JL + . For our choice of variables = {a, 

p^ = {—p a ,P<t>} metric has the form (in the units 4ttG/3 = 1) g^ v 



/-A \ 

n 



„3 



SO 



that g = a 4 , p a = i (JA + Aj and the Hamiltonian is 

~ 1 / 2 1 1 a2 \ pi 1 d d Id 2 

Explicit expression for the wave function satisfying Hqi/j = is 
il> k (a,<f>) = a ±m c ikrt> . 

Exactly as in the case of the Klein-Gordon equation we should choose only the 
positive frequency solutions j^]. Such wave function corresponds to the definite choice 
of the boundary condition for the minisuperspace: the wave function is formed by the 
modes bounded on k and only ingoing in a singularity. 

Thus the wave packet 

— i\k\ 

V»(o,0) = / c(k) a e ikrp dk (13) 
3 y/4ir\k\ 

will be normalized by 

(dip dip* 



HI 



I v£ r ~ it*) d(p= I c *^)< k ) dk = l ( i4 ) 

The proper time evolution of the operators 
results in 

3 3 

p^r) = 0, p^(r) = const, (a 3 (r))' = - (p a a + ap a ) , ~(p a a + ap a )' = -9H . (15) 

Set of equations for the Heisenberg operators and the constraint equation for the 
states would be considered as a tool to evaluate the mean values of observables. The 
obstacle is that the operators are hermite relatively an integration over ^gd n x, while 
due to constraint the states can not be normalized in such a way and are normalized in 
the Klein-Gordon style. 

Our recept consists in enforcing the constraints on the equation of motion for the 
Heisenberg operators at r = 0. First let us remind the Dirac quantization procedure 
PU] and return to the classical picture for this goal. According to Dirac, besides the 

2 p2 _____ 

primary constraint $i = — ^f + ^f (see (jTUjl ). we have to set some additional gauge fixing 
(secondary) constraint, which can be chosen in our case as $2 — a ~ const, because the 
hyperplane a = const is chosen earlier for the normalization of the wave function in 
the Klein-Gordon style. In contrast to the usual formalism El I2U 122]) we impose 
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constraints $1 = and $2 — 0, which have to be satisfied only at hyperplane r = 0. 
Besides the ordinary Poisson brackets 

B 3A_dB_ _ dA_dB_ 
dp^ dxt 1 dx^ dp^ ' 
the Dirac brackets are introduced 

{A, B} D = {A, B} - {A, ^}(C-%{^, B}, (17) 

where C is the nonsingular matrix with the elements CV,- = and is the 

inverse matrix. Quantization consists in the substitution of the commutators in the 
Dirac brackets and the replacement of the variables by the operators: 

[fi,f)'] = -i{r,,r,'} D . (18) 

Here r) implies the set of the canonical variables p^-,x v . 

The Heisenberg equations of motion have to satisfy the constraints at initial moment 
t = and the operators obey the commutation relations obtained from the Dirac 
quantization procedure. Direct evaluation gives 

[p«(0),a(0)] = 0, [p (O),a(O)] = 0, 

[fr(0), 0(0)] = -i, [p o (0), 0(0)] = -i - (19) 

^(0)^(0) 

We have to solve the equations (fTo"j) with the given initial commutation relations. 
This changes the canonical commutation relations between the Heisenberg operators at 
r / 0. Therefore we shall call theirs as the quasi- Heisenberg operators. In fact we 
appeal to the structure of the classical theory and re-quantize the equations of motion. 

One can see that the commutation relations ("T"*"|) can be satisfied through 

a(0) = const = a, p^O) = p^, p a (0) = \p^\/a, 0(0) = 0, (20) 

where p^ = Variable a is c-number now because it commutes with all operators 

j. Solutions of the equations (fT5")) are 

P4>i. T ) = Ph <K r ) = <P+ 777~~"~T m (« 3 + 3 l^l r ) - In a, 



Pa{T) = (a^ + 3\l\ry^ a3(T) =a3 + 3| ^ |T ' 

(21) 

We imply that these quasi-Heisenberg operators act in the Hilbert space with the 
Klein-Gordon scalar product. 

Mean value of some operator A(t) can be defined as [2*3] : 

<9-0(a, 0) 



< A(t) >= ia J (ifi*{a, (f))A(r, (f>,p,p, a 
dip* (a, 0) 



da 



da 

A(t, 0,p^,a)^(a, 0) )d<p 



(22) 

a^0 
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However we shall use another definition: 

i d 



< A(t) >=iaj U\a,4>)D l ±A{T)D-*—il){a,<j)) 

(23) 



-?-iP*(a,(f>)D-U(T)D^(a, 
oa / a ^o 

where operator D = — ^ + 2a 6 V((p) (since a — ► the l^-term can be omitted in the 
expression for D). The mean value (J23|) is particular case of that suggested in Ref. [24J, 
where an one-particle picture of the Klein-Gordon equation in the Foldy-Wouthausen 
representation is considered. The advantage of this normalization can been seen in the 
momentum representation of <p variable, where p^ = k and <p = iJ^. Eq. (|2*2*|) gives 

< A(t) >=\J a iW ^\k\c*{k)A{r^,k,a)a^^ldk 



(24) 



/ a 4|fc| ^^A(r,0,fc,a)a- J|fe| v /|fc|c(fc)^ 
2 wife' 



while Eq. (|23j) leads to 

< A{t) >= J a i l fe l C *(A;)i(r,0,fc,a)a- J l fc lc(A;)rfA; 

which is similar to the ordinary quantum mechanics and certainly posses hermicity. 
Evaluation of the mean value over the wave packet gives 

< a 3 (r) >= 3r / \k\\c(k)\ 2 dk. (25) 



If we assume that there is no external "physical time" in Universe (proper time r hardly 
can be considered as the measurable observable), then we must take some quantity, for 
instance, < a 3 (r) > as the "physical time". 

Next interesting quantity is the mean value of the scalar field < <f)(r) >: 

<0(r) > a = / a l ^c*(k) (i-^- + -^-\n(a 3 + 3\k\r] - In a ) a^cf^dk 
J V ok 3\k\ v ' \k\ J 

= /( c *( A; )^ c ( A; ) + ^ m ( a3 + 3 l A; l r )l c ( A; )| 2 )^- ( 26 ) 

Brackets < ... > a with the index a in (J26|) mean that we do not set a equal to 
zero yet (compare with Eq. (|24j)). A remarkable property of Eq. ()26|) is that the term 
— j|| In a cancels the term arising from the differentiation: a*' fc 'z J^a~*' fc ' = ||ylna. Thus 
we may get to a — > and obtain 

<kr) >= J ^ln(3|fc|r)|c(A;)| 2 + c*(^c(^^. (27) 

Cancellation of the terms divergent under a — > in the mean values of the quasi- 
Heisenberg operators is a general feature of the theory, and give us possibility to evaluate, 
for instance, 

<0 2 (r)>= |Qln 2 (3|A;|r)|c(A;)| 2 -c*(A;)^c(fc))^. (28) 
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One should not confuse the divergence at a — > arising under evaluation of the 
mean values with the singularity at r — > 0. The mean values of operators, which are 
singular at r — > in classical theory remain singular also in the quantum case. The way 
to avoid a singularity is to guess, that Universe was burn not from a point but from a 
" seed" of ao-" size" . Then in the expression for mean value we have to assume a — > 
instead of a — > 0. This puts a question about underling theory giving size of the seed. 

One more kind of the infinity can be found in Eqs. ([27)1 . ([28)1 . For c(k), which 
does not tend to zero at small k the mean values of 0(r) and 2 (t) diverge. This is a 
manifestation of the well known infrared divergency of the scalar field minimally coupled 
with gravity. Thus not all possible c(k) are suitable for the construction of the wave 
packets. 

Let us consider Hamiltonian, containing the cosmological constant Vq. 

H = H + a 3 V . (29) 
Explicit solution for the wave function Hip = has the form 

18^ W1 i|*| w ,vTOs^<* 



MaA) = [yj r(i - -y)^-f(Y a )e ' (30) 

where T(z) is the Gamma function and J^z) is the Bessel function. The wave function 
(1371 tends to tH^e^l + 0(a 6 )) asymptotically under a 0. If we assume that the 
operators of all physical quantities are local on the a variable or can be approximated by 
the local ones, then for an evaluation of the mean values we may always build the wave 
packet from the functions a _l ' fc 'e ifc<? \ This supply ideology of the usual nonrelativistic 
Heisenberg picture which says that all dynamics is contained in the Heisenberg operators 
and only knowledge of the wave function is required. The argumentation holds for any 
potential V(4>), because it contributes into the Hamiltonian as a term multiplied by a 3 . 
Equations of motion obtained from the Hamiltonian (|29|) are 



3 

(a 3 (r))" = -{p a a + ap a ); (31) 

^{p a a + ap a y = 9V a 3 - 9H ; (32) 

9(V a 3 - H y = 18V ^(p a a + a Pa ). (33) 
The additional term a 3 Vo does not change relations ()19j) required for the re-quantization 



procedure. Only expression for p a (0) changes in (f2TIj) : p a (0) = y + 2VoCt 4 
Finally we arrive to 



a\r) = a 3 + 3|^|^^1 + a * (cosh(r ^8V~ ) - 1). (34) 

V to Vq 

Evaluation of the mean values according to Eq. (|2*3*|) leads to 



< aHr) >= 5in ' l( ^f T> / 

fsinh(3 v / 2l^r) 

< a\r) >= A } > I lr\r(k)\^k. 
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This shows, that the dispersion ^ <a <? > a 3>'^ ~ does not depend on r, exactly as in the case 
of the free relativistic particle. Thus in this model the evolution of Universe remains 
quantum during all time. This results from the absence of any scale length (but not 
from an ambiguity of the wave function normalization; see, for instance, Ref. J25J). 
Such a length appears if we take V{(j)) = One may suggest that the expansion of 

Universe remains quantum until the moment, when the scaling factor approaches the 
Compton wave length 1/m. When < a(r) > becomes greater than 1/m the expansion 
can be described classically. To see this explicitly we have to find the quasi-Heisenberg 
solutions with the above potential. 



4. Operator equations for the quadratic inflationary potential 

As it has been discussed the quantization procedure is reduced to the quantization of 
the equations of motion i.e. considering them as the operator equations, which have 
to be solved with the initial conditions satisfying to the constraint at r = 0. For the 
Hamiltonian 

we have the equations 

3 • o (i 3 o 
a = — am 1 — aqm , 

= -3-0 - gdj (36) 
a 

and the constraint 

- a 2 a + 0V + a 3 g<j) 2 = 0. (37) 

The point means the differentiation over r. After quantization Eqs. (|36|) lead to the 
equations for the quasi-Heisenberg operators, which have to be solved with the operator 
initial conditions 

9 i, n \ P*(°) 1 ( ■ 9 \ 



a(0) = a, 0(0) = 0, p^(0) = -i— , 0(0) 



a 3 (0) a 3 



—i- 



a(O) = a(OW0 2 (O)+£0 2 (O) 



\ 



If . d N 2 



a 



l 



i— + ga 2 (f) 2 . (38) 



According to our ideology the operator constraint (|S7|) is satisfied only at r = 0. The 
ordinary problem of the operator ordering arises, because in the general case the quasi- 
Heisenberg operators are noncommutative. The problem becomes more transparent if 
we change the variable a = In a: 

a + -a - -g(j) + -<f> =0, 

+ - (&<f) + (f>aj +gdl = 0. (39) 
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Figure 2. Mean values of the square of the scalar field and the logarithm of the 
Universe scale factor calculated in the zero (1) and first orders (2) on g = 1, for the 
wave packet c(fc) = -4| (^) 1 ^ 4 fc 2 exp(— fc 2 ) - Fig. (a), (b), and for the wave packet 

c(k) — (f) 1 ^ 4 exp(— k 2 — + 2) - Fig. (c), (d). Each of the triplets of the curves 
in Fig. (b), (d) corresponds to the < d(r) > —y/D(r) (lower), < &{t) > (medial), 
< d(r) > + v /D(t) (upper), where D =< a 2 (r) > - < a(r) > 2 . 



This system has to be solved with the initial conditions 0(0) = <p, a(0) = Ina, 

0(0) = ^3 ( — ^(0) = y^^ 2 + ^ ( — • We have used the symmetric ordering in 
Eq. flUD- 

The operator equations under consideration can be solved within the framework 
of the perturbation theory in the first order on interaction constant. The solution 
in analytical form is given in Appendix. The calculations of the mean values based 
on (j24j) have been carried out with the two kinds of the wave packets: c(k) = 

^ k 2 exp(— k 2 ) and c{k) = e 2 exp(— k 2 — -p), which were normalized as 

\c{k)\ 2 dk = 1. Because of the infrared divergency we must take only wave packets 
with vanishing c(k) at k — > 0. The results of the calculations are shown in Figs. [21 
The packets under consideration are symmetric thus the mean value of the scalar field 
vanishes. In the first order on g we can consider only an early beginning of the inflation. 
At this stage the square of the the scalar field grows due to initial "kinetic energy" . The 
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logarithm of the scale factor grows linearly. 

As it has been discussed in the previous section in the zero-order on g, the relative 
dispersion of a 3 (r) does not depend on r so that the dispersion of d(r) = In a(r) is 
constant. In the first order on g we might expect that the dispersion of the scale factor 
should decrease and the quantum universe will approach to the classical one. However, 
we see an opposite picture: the dispersion of a grows with r. However, as it is suggested 
in Ref. |2Z1 some mechanism of the classical world appearance has to exist. 



5. Wigner-Weyl evolution of the minisuperspace 



The analysis of the inflation at its late stages requires a numerical consideration of 
the operator equations that can be realized within the framework of the Weyl-Wigner 
phase-space formalism [23J. Let us remind that in this formalism every operator acting 
on if variable have the Weyl symbol: W[A] = A(k,(j>). The simplest Weyl symbols 
reads: W[— i-jn] = k, W[0] = 0. Weyl symbol of the symmetrized product of operators 
has the form 



W[-(AB + BA)} = cos - 



h d d h d d 



A{k u <j> x )B{k 2 



k\=k2=k, ip 1= (f)2=(p 



hdk 2 2d(j) 2 dk 1 

where the Plank constant is restored only to point up to which order we will expand 
the cosine in the next. This has no direct physical meaning because the effects of the 
quantum mechanics are contained as in the Weyl symbols, so in the Wigner function, 
and the last one has no limit at h — > in the general case f2*3l . 

Let us consider the Weyl transformation of Eqs. (jHEJ) and expand the Weyl 
symmetrized product of operators up to second-order in H. This results in: 

d 2 a + \([d T af + ^(d k d,d T af - ^(^ T a)(^ r a)) +\{[d rV f + ^(d k d,d T pf 

-^(<^)(^)) - \ g (<? + f (d k d, v ) 2 - ^)(^)) = o, 



h 2 



n 2 



d 2 (p + 3[d T ad T (p + ^(9 fc ^9 T a)(9 fc ^9 r v9) - '^-(d 2 k d T a)(dld T (p) 

-%{dld T a){dldM 



9¥ 



0. 



(40) 



where a(k, 0, r) and ip(k,(p,r) are the Weyl symbols of the operators d(r) and <p(r) 
respectively. These equations have to be solved with the initial conditions at r = 0: 



a (k, 0, 0) = In (a) , d T a (k, 0, 0) = W 



1 d 2 



+ 



if (k, T, 0) = 0, d T ip (k, T, 0) 



k 



(41) 



Weyl symbol of the square root 26 can be expressed as 
W go 1 



a 6 <90 J 



1/4 



9 

7r 1 / 2 a 3 / 2 Jo 



t 



-1/2 



cxp 



ga 



6^2 



+ k 2 



V9 ( 



tanh(t) 
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xsech(t) [ ga sech(t) 2 + tanh(t) ) dt, 



(42) 



v Vd a3 / 
but since the mean values are calculated in the limit a — > 0, we can take simply 
d T a (k, (f>, 0) = ^ (this will not change the mean values). 

State of the Universe in this approach is described by the Wigner function p(k, <fi), 
which is constructed on the basis of definition (J23j) . That gives: 

3 dip (0 - f ) ' 



p (k, 0) = ia 

d 2 J9f (0 + l) 



9 2 



^(0+2 



a 2 



5 2 



9a 



da 
e iku du, 



e iku du 



(43) 



or in the momentum representation of the wave function corresponding to Eq. (|13jl: 

p(Jfe, </>) = -[ c*(2k - g)c(g)a- i|9l+i|2fc - 9l e 2i(9 - fc) ^g. (44) 



7T 

As a result of a — > 0, both Weyl symbols and Wigner function diverge. In particular, 
when a — > the Wigner function becomes more and more oscillating. However the 
divergences cancel each other in the expectation values constructed in the ordinary way. 
For instance, expectation values of a and its square are: 



(a(r)> 



a 



2 V) 



dkd(pa(k, <p, r)p (k, <f))\ a ~,o 

h 2 



(dkd^a) - — (dj,a)(d k a) p{k,4>) 



J dkdcf) yy 2 + ^ v _ n ,„^„ / ^ 

An example of the Wigner function for the wave package with c (k) = 

771 (f) 4 /c 2 exp (—k 2 ) is shown in Fig. El 

As a result of numerical solution of Eqs. ([40)1 . we have obtained the evolution of the 
operators expectation values and their dispersions. Some results are shown in Figs. (0J- 
EJ). We have considered three wave packages providing suppression of the infrared (and 
certainly ultraviolet) divergence. For our rather illustrative calculations, the "cuts off" 
of the wave packets are chosen of order unity i.e. at the Plankian level. In principal, the 
"cuts off" have to be consistent with the sub-Planskian physics eliminating divergences 
at fundamental level. In other words we describe Planskian physis here, but all sub- 
Planskian physics should be contained in the functions of c(k). 

It is of interest to compare the evolution numerically obtained for the arbitrary 
coupling parameter g with that obtained from the analytical operator evaluation in the 
first-order expansion on g (see previous Section). We consider two different cases for 
Eqs. (JIQJ): i) h — (quasi-classical evolution) and ii) % = 1 (quantum corrections; note, 
that in the previous sections we didn't use any expansion on %). 

In addition to results of previous Section, we can see the deceleration of the 
initial exponential expansion of the universe and the transition to the post-inflationary 
scenario. In the quasi-classical case the evolution of expectations can be considered as 
a good approximation to the classical scenario (curves 1 and dashed (or dotted) curves, 
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-2 2 

k 




Figure 4. Evolution of (a) , (<p 2 ), dispersion a (a) — y (a 2 ) — (a) 2 (black curves) 

and relative dispersion a (a) / (a) (gray curves) for c(k) — -j= (§) 4 fc 2 exp(—fc 2 ). 
Ti = (1), 1 (2). For comparison we give the classical evolution with the same initial 
conditions (initial k is (dashed) and y (k 2 ) (dotted), initial 4> =\J (^ 2 ))- a = 10 -4 . 



respectively). However, the inflaton oscillations are smeared in the quasi-classical case 
due to presence of the modes with different k. Choice of the more "cold" wave package, 
i.e. the package with the suppressed high-frequency modes, enhances the inflation (Fig. 
EJ). Some enhancement of inflation results also from the removing of the low-frequency 
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<a> <$ 2 > a(a) a(a)/<a> 




1 10 20 30 1 10 20 30 1 10 20 30 



Figure 6. As in Fig. Hbut for c(k) = e 2 (f)* exp (-fc 2 - 1/fc 2 ). 

components from the wave package (Fig. 

The quantum corrections (second-order in Ti) essentially affect the evolution at the 
initial stage. As it was demonstrated in previous Section, there is the stage when the 
scalar field rolls away the minimum of potential due to kinetic energy of the wave package 
(some analog exists also in the classical case, see dotted curves in Fig. EJ). As a result, 
the inflation intensifies. However, in any case the inflation comes to the end. 

It is important, that in all cases the relative dispersion reaches some maximum and 
then decreases. However, the relative dispersion does not vanish but approaches some 
asymptotical value, which does not differ essentially for the different wave packages. It 
is astonishing that the asymptotical relative dispersion is similar in quasi-classical and 
in quantum-corrected cases. This suggests that some mechanism of the classical world 
appearance is required. Such mechanism is absent in our simplest model. We surmise 
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that this can be some decoherence produced by the additional degrees of freedom. 
6. Conclusion 

Quantum evolution of the Universe originated from the some fluctuation of the scalar 
field (wave packet) has been considered. 

Re-quantization procedure for the Heisenberg operators has been introduced to 
compensate the loss of hermicity arising from the unboundedness of the Universe wave 
function along a variable. 

Mean values of the operators corresponding to the observables, which are singular 
at t —>■ on a classical level, remains singular also in a quantum case. 

If the scale like the mass of the scalar field is not introduced, the evolution of 
Universe remains quantum eternally in the sense that the relative dispersion of the scale 
factor does not decrease with the increasing r. 

We have solved the equations for the quasi-Heisenberg operators in the first order 
on g for the inflation potential V(4>) = g^r. In the first order on g the dispersion of the 
logarithm of the scale factor increases with the increasing r. 

Numerical calculations on the basis of the Weyl-Wigner phase-space formalism have 
demonstrated the exit from the inflation with the approaching of the relative dispersion 
to some constant but essentially non-zero value. This value does not depend noticeably 
on the universe's wave package. 

We are forced to establish that we fail to come to "classical" Universe (i.e. the 
Universe with small dispersions) from the quantum state possing large dispersions. 
Certainly we can take the state possessing small dispersion and thus obtain the 
" classical" Universe, but it is not so interesting. More interesting is to find some peculiar 
mechanism of the classical world appearance. 

Appendix A. Solution of the equations for the quasi-Heisenberg operators 
in the first order on g 

Here we listed the solution of the operator equations for the quasi-Heisenberg operators 
in the first order on g. Evaluation is done in the momentum representation of the 
variable, where the initial conditions have the form: a(0) = In a, 0(0) = -4-, 

0(0) = = (ijQ, d(0) = V30 2 + g = \jg (^) 2 + 5- In the first order on 9 the 
square root can be extracted giving 

The solutions of have been obtained in the framework of the perturbation theory 
using a computer algebra and have the form: 

x 7 k f a 3 + 3\k\r \ f F 2 2 + 20F 2 + 2 F 2 ^0 + 0^ \ 

^ T) = ^ + m \ ^~) +9 [ 4 + — 2 — +Fo ) 



18 



= 3 ln ( a + 3|fc|r) + # I + + A) I , (A.2) 



where functions F 2 , D 2 , F ly . . . are 

3a 6 VFr 2 -6a 3 fc 2 t 3 - 9 (k 2 ) 1 t 4 



F 2 (k,T,a) = 



-4a 6 k + 36 fc 3 r 2 



Fl (fc ' a) = -108^^ + 972 fc5 T 2 (~6 a 9 fc VP r + 27 a 6 fc 3 r 2 - 18 a 3 fc 3 ^r 3 - 27 fc 5 r 4 
-2 a 6 k (a 6 - 9 fc 2 r 2 ) ln(a 6 ) + 2 a 6 VP (-a 6 + 9 k 2 r 2 ) ln(l - ™±) + 2 a 12 VP ln(l + ™±) 
-18 a 6 (fc 2 ) 1 r 2 ln(l + "*Al) - 2 a 12 fc ln(l + + 12 a 9 fc VP r ln(l + 

+36 a 6 fc 3 t 2 ln(l + - 108 a 3 fc 3 #r 3 ln(l + '^^-) - 162 fc 5 r 4 ln(l + 

+2 a 12 fc ln(a 6 - 9 fc 2 t 2 ) - 18 a 6 fc 3 r 2 ln(a 6 - 9 fc 2 t 2 ) 

F (fc, T, a) - — — — — L ( _ 6 a 9 fc 2 T _ gl fl 6 /J.2N I T 2 + 19g fl 3 fc 4 T 3 + 397 fc 4 ^ ^ + /_ g fl 12 ^ 

1944 fc" 5 (— a D + 9 Ar t a ) \ 

+72 a 6 (fc 2 ) 1 r 2 ) ln(a 6 ) + (-8 a 12 fc + 72 a 6 fc 3 r 2 ) ln(l - ^Al) + (8 a 12 fc - 72 a 6 fc 3 r 2 ) ln(l + *?Al) + 

(-14 a 12 Vfc 2 " + 12 a 9 fc 2 r + 144 a 6 (fc 2 ) 1 r 2 - 108 a 3 fc 4 r 3 - 162 fc 4 #t 4 J ln(l + - 6 (-a 6 + 9 fc 2 r 2 ) 

x (a 6 Vfc 2 " + 6 a 3 fc 2 r + 9 (fc 2 ) 1 r 2 ) ln(l + + 8 a 12 VP ln(a 6 - 9 fc 2 r 2 ) - 72 a 6 (fc 2 ) 1 r 2 ln(a 6 - 9 fc 2 r 2 )) 

a 3 / „ r-r loI , a 12 fc 2 a 12 kV¥ 3a 9 fc 2 3a 9 fcVF „ B sViF l2 

+ 7T^-r^ - 2 a 3 fc Vfc 2 - 18 fc 3 r 7 + t + ~ ? - 6 a 6 s~ fc 2 

216fc 6 V (a 3 -3fcr) 3 (a 3 -3fcr) 3 (a 3 -3fcr) 2 (a 3 -3fcr) 2 

x( a 3 -3fcr)- 1+ ^+3 a 6 -^fc^( a 3 -3fcr)- 1+ ^ + 18a 3 -^fc 3 r( fl 3 -3fcr)- 1+ ^ ! + k - 3 

V ' 17 V ' (« 3 + 3fcr) 3 

a^fcVfc 2 3a 9 fc 2 S^fcVfc 2 e^+^fc 2 S^+^fcA/fc 2 18 a fc 3 r 

+ 7= + 7=" + 



(a 3 + 3fcrr (a 3 + 3fcrr (« 3 + 3fcrr (a 3 + 3fcr) *±# (a 3 + 3 fc T )^ ( a 3 + 3 fcr)^ 
+ 9 ° " ^ + ^=I-9(fc 2 ) f T^l-^Al^ " -18a 3 fcVF ln(a) + 3a 3 fc (-fc + Vfc 2 ") ln(a 3 -3fcr) 



(a 3 + 3fcT)^ 
+3a 3 fc 2 ln(a 3 + 3fcr) + 3a 3 kVk? ln(a 3 + 3fcr) 



D 2 {k,T, a) 



a 



fc 2 r - 6 a 3 (fc 2 ) 3 r 3 - 9 fc 4 r 4 
2a 6 fc 2 - 18fc 4 r 2 



Dl {k > T ' a) = 972 fc 3 (- fl 6 + 9fc 2 r 2 ) ( 108 ° 9 ^ T + 486 ^ 1 T ' " 1620 ° 3 ^ T ' " 2430 ^ ^ r ' 

-36 a 12 ^ ln(l + _ 2 16 a 9 fc 2 r ln(l + + i 944a 3 fc 4 r 3 ln(l + '^^-) 

a 6 a A a A 

+2916 fc 4 Vk^r 4 ln(l + ^f 1 )) ; 
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D (k,T,a) 



972 fc 3 (-a 6 + 9fc 2 r 2 



-30 a 9 kVk^r - 81 a 6 fc 3 r 2 + 342 a 3 k 3 Vk? t 3 + 513 fc 5 t 4 



40 a 12 k ln(l + ^ll) + 60 a 9 fc Vfc 2 r ln(l + 1^1) _ 540a 3 fe 3 #r 3 ln(l + 



— — 2 

-•->,„,, , 3v ^)-6a 12 fcln(l + ^-^) -36a 9 fcVPrln(l 



-810fc 5 r 4 ln(l 



324 a 3 k 3 Vk^r 3 ln(l + + 4 86 fc 5 r 4 ln(l + 



216 fc 6 



3Vfc^T 

a 3 } + 
- 8fc 2 



a 9 fc 2 



a 9 k Vk? 



3a 6 k 2 3a 6 kV¥ 



., , 2 +— : 2 + 6a 3 ~^~ fc 2 (a 3 - 3fcr) 



-i+- 



(a 3 -3fcr) J (a 3 -3fcr) J (a 3 -3fcrr (a 3 -3fcr) 



18fc 3 r (a 3 -3kr) 



3 Vk^ 



9(fc 2 ) i r( Q 3 -3fcr)' 1 ' 



a 9 /fc 2 



a 9 fc V^fc 2 



(a 3 + 3fcr) 3 (a 3 + 3fcr) 3 



a k a k 

3a 6 fc 2 3a 6 kV¥ — " fc 2 3a^ — "JfeVF 18a^fc 3 r 



(a 3 + 3fcr) (a 3 + 3fcr) (a 3 + 3A:r) 



9 — 3 fc^ri-^ 



(a 3 + 3fcr) fc 

„2 , / 3 



(a 3 + 3fcT) — ~ (a 3 + 3fcr)" 
18 fc 2 ln(a) + 3 fc (fc - ln(a 3 - 3 fc r) 



The output has been generated from MATHEMATICAL \f}£ means \k\ 
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